all the pixels along the spectral-disper- 
sion dimension would be summed to 
obtain the value of the cross-correlation 
(plus background) . 

Such on-chip cross-correlation could 
be performed rapidly because the ana- 
lytical function could be statically pro- 
grammed into the APS array and the 
multiplications could be done simulta- 
neously or nearly so. All of the addi- 
tions could be done simultaneously by 
means of a single binning instruction. 
The charge wells of all the pixels could 
be connected simultaneously, collect- 
ing all the charge outputs from multi- 
plication operations into one “super- 
pixel,” the single readout value of 
which would constitute the cross-corre- 
lation value for the given analytical 
function. For an instrument in which 
the APS rows were aligned along the 
spectral-dispersion dimension and in 
which the image of a spectrograph slit 
was aligned along the pixel columns 
and spanned multiple pixel rows, it 


would be possible to perform simulta- 
neous cross-correlations for multiple 
target species by applying, to each pixel 
row, the analytical function correspond- 
ing to one of the target species. A sepa- 
rate readout would be needed for each 
target species. 

In the other hardware implementa- 
tion, cross-correlations would be com- 
puted externally to the APS array. The 
multiplications and additions would be 
performed in pipeline fashion. If the 
APS-array outputs were analog, then 
programmable analog signals repre- 
senting the analytical functions would 
be synthesized in phase with the corre- 
sponding stream of analog APS-array 
outputs and the multiplications and ad- 
ditions would be performed by rela- 
tively inexpensive, commercially avail- 
able analog mixing and filtering 
circuits, respectively. If the APS-array 
outputs were digital, the cross-correla- 
tions could be computed by a digital 
signal processor. Ordinarily, the analog 


approach would be preferable because 
the analog operations can be per- 
formed much more rapidly than can 
the corresponding digital multiplica- 
tions and additions. 

This work was done by Gregory Bearman, 
Michael Pelletier, and Suresh Seshadri of 
Caltech for NASA’s Jet Propulsion Labo- 
ratory. Further information is contained in 
a TSP (see page 1 ). 

In accordance with Public Law 96-517, 
the contractor has elected to retain title to this 
invention. Inquiries concerning rights for its 
commercial use should be addressed to: 
Innovative Technology Assets Management 

JPL 

Mail Stop 202-233 
4800 Oak Grove Drive 
Pasadena, CA 91109-8099 
(818) 354-2240 
E-mail: iaoffice@jpl.nasa.gov 
Refer to NPO-30912, volume and num- 
ber of this NASA Tech Briefs issue, and 
the page number. 


Prioritizing Scientific Data for Transmission 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


A software system has been developed 
for prioritizing newly acquired geologi- 
cal data onboard a planetary rover. The 
system has been designed to enable effi- 
cient use of limited communication re- 
sources by transmitting the data likely to 
have the most scientific value. This soft- 
ware operates onboard a rover by ana- 
lyzing collected data, identifying poten- 
tial scientific targets, and then using 
that information to prioritize data for 
transmission to Earth. Currently, the sys- 
tem is focused on the analysis of ac- 
quired images, although the general 
techniques are applicable to a wide 
range of data modalities. Image prioriti- 
zation is performed using two main 


steps. In the first step, the software de- 
tects features of interest from each 
image. In its current application, the 
system is focused on visual properties of 
rocks. Thus, rocks are located in each 
image and rock properties, such as 
shape, texture, and albedo, are ex- 
tracted from the identified rocks. In the 
second step, the features extracted from 
a group of images are used to prioritize 
the images using three different meth- 
ods: (1) identification of key target sig- 
nature (finding specific rock features 
the scientist has identified as impor- 
tant), (2) novelty detection (finding 
rocks we haven’t seen before), and (3) 
representative rock sampling (finding 


the most average sample of each rock 
type). These methods use techniques 
such as K-means unsupervised cluster- 
ing and a discrimination-based kernel 
classifier to rank images based on their 
interest level. 

This program was written by Rebecca Cas- 
tano, Robert Anderson, TaraEstlin, Dennis De- 
Coste, Daniel Gaines, Dominic Mazzoni, Forest 
Fisher, and Michele Judd of Caltech for NASA’s 
Jet Propulsion Laboratory. Further informa- 
tion is contained in a TSP (seepage 1 ). 

This software is available for commercial li- 
censing. Please contact Don Hart of the Cali- 
fornia Institute of Technology at (818) 393- 
3425. Refer to NPO-40265. 


Determining Sizes of Particles in a Flow From DPIV Data 

The same equipment would be used to measure sizes as well as velocities. 

John H. Glenn Research Center, Cleveland, Ohio 


A proposed method of measuring the 
size of particles entrained in a flow of a 
liquid or gas would involve utilization of 
data from digital particle-image ve- 
locimetry (DPIV) of the flow. That is to 
say, with proper design and operation of 


a DPIV system, the DPIV data could be 
processed according to the proposed 
method to obtain particle sizes in addi- 
tion to particle velocities. As an addi- 
tional benefit, one could then compute 
the mass flux of the entrained particles 


from the particle sizes and velocities. 

As in DPIV as practiced heretofore, a 
pulsed laser beam would be formed into 
a thin sheet to illuminate a plane of in- 
terest in a flow field and the illuminated 
plane would be observed by means of a 
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Read in 12-bit image file. 

] 

Apply global threshold 
to eliminate image noise. 

I 

Scan image for nonzero 
pixels and group nonzero 
neighbors into images of 
particles. 

\ 

Determine glare-spot 
positions and calculate 
particle diameters. 

\ 

Calculate positions of 
particles, the distribution of 
particle sizes, and the 
mean and standard 
deviation of particle sizes. 


The Algorithm of the Proposed Method would 
process DPIV data to extract information on par- 
ticle sizes. The size information would be in ad- 
dition to the velocity information extracted 
from the data by use of previously developed 
DPIV software. 

charge-coupled device (CCD) camera 
aimed along a line perpendicular to the 
illuminated plane. Unlike in DPIV as 
practiced heretofore, care would be 
taken to polarize the laser beam so that 
its electric field would lie in the illumi- 
nated plane, for the reason explained in 
the next paragraph. 


The proposed method applies, more 
specifically, to transparent or semitrans- 
parent spherical particles that have an 
index of refraction different from that of 
the fluid in which they are entrained. 
The method is based on the established 
Mie theory, which describes the scatter- 
ing of light by diffraction, refraction, and 
specular reflection of light by such parti- 
cles. In the case of a particle illuminated 
by polarized light and observed in the 
arrangement described in the preceding 
paragraph, the Mie theory shows that the 
image of the particle on the focal plane 
of the CCD camera includes two glare 
spots: one attributable to light reflected 
toward the camera and one attributable 
to light refracted toward the camera. 
The distance between the glare spots is a 
known function of the size of the parti- 
cle, the indices of refraction of the parti- 
cle material, and design parameters of 
the camera optics. Hence, the size of a 
particle can be determined from the dis- 
tance between the glare spots. 

The proposed method would be im- 
plemented in an algorithm that would 
automatically identify, and measure the 
distance between, the glare spots for 
each particle for which a suitable image 
has been captured in a DPIV image 
frame. The algorithm (see figure) 
would begin with thresholding of data 
from the entire image frame to reduce 
noise, thereby facilitating discrimina- 
tion of particle images from the back- 
ground and aiding in the separation of 
overlapping particles. It is important 
not to pick a threshold level so high 
that the light intensity between a given 
pair of glare spots does not fall below 
the threshold value, leaving the glare 
spots disconnected. 

The image would then be scanned in a 
sequence of rows and columns of pixels 
to identify groups of adjacent pixels that 
contain nonzero brightnesses and that 


are surrounded by pixels of zero bright- 
ness. Each such group would be assumed 
to constitute the image of one particle. 
Each such group would be further ana- 
lyzed to determine whether the image 
was saturated; saturated particle images 
must be rejected because the locations of 
glare spots in saturated images cannot 
accurately be determined. Within each 
unsaturated particle image, the cen- 
troids (deemed to be the locations) of 
the glare spots would be determined by 
means of gradients of brightness distrib- 
utions and three-point horizontal and 
three-point vertical Gaussian estimates 
based on the brightness values of the 
brightest pixels and the pixels adjacent 
to them. If the brightness of a given par- 
ticle image contained only one peak, 
then it would be assumed that a second 
glare spot did not exist and that image 
would be rejected. 

Once the centroids had been esti- 
mated for all particle images for which 
it was possible to do so, the positions of 
the particles and the distances between 
their centroids would be computed. As 
described above, the size of each parti- 
cle would then be computed from the 
distance between its centroids. Finally, 
the distribution, mean, and standard 
deviation of sizes would be computed 
for the collection of particle images 
that survived to the final stage of the 
centroid-estimation process. 

This work was done by M. P. Wernet of 
Glenn Research Center and A. Mielke 
and J. R. Kadambi of Case Western Reserve 
University. Further information is con- 
tained in a TSP (see page 1 ). 

Inquiries concerning rights for the com- 
mercial use of this invention should be ad- 
dressed to NASA Glenn Research Center, 
Commercial Technology Office, Attn: Steve 
Fedor, Mail Stop 4-8, 21000 Brookpark 
Road, Cleveland, Ohio 44135. Refer to 
LEW-17340. 


Faster Processing for Inverting GPS Occultation Data 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


A document outlines a computa- 
tional method that can be incorpo- 
rated into two prior methods used to 
invert Global Positioning System 
(GPS) occultation data [signal data ac- 
quired by a low-Earth-orbiting satellite 
as either this or the GPS satellite rises 
above or falls below the horizon] to 
obtain information on altitude-depen- 


dent properties of the atmosphere. 
The two prior inversion methods, 
known as back propagation and canon- 
ical transform, are computationally ex- 
pensive because for each occultation, 
they involve numerical evaluation of a 
large number of diffraction-like spatial 
integrals. The present method involves 
an angular-spectrum-based phase-ex- 


trapolation approximation in which 
each data point is associated with a 
plane-wave component that propa- 
gates in a unique direction from the 
orbit of the receiving satellite to inter- 
sect a straight line tangent to the orbit 
at a nearby point. This approximation 
enables the use of fast Fourier trans- 
forms (FFTs) , which apply only to data 
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